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Abstract 

Adsorption at a 1- dimensional planar substrate equipped with a 
localized chemical inhomogeneity is studied within the framework of a 
continuum interfacial model from the point of view of interfacial mor- 
phology and correlation function properties. Exact expressions for the 
one-point and two-point probability distribution functions Pr(lr) and 
Pr 1 ,r 2 (^r 1 , £r 2 )> being the interface position above a fixed point T of 
the substrate, are derived for temperature corresponding to the inho- 
mogeneity's wetting transition. It is demonstrated that in the limit of 
macroscopic inhomogeneity's size the net effect of the remaining ho- 
mogeneous parts of the substrate on the interfacial morphology above 
the inhomogeneity is exactly equivalent to appropriate pinning of the 
interface at its boundaries. The structure of the average interfacial 
morphology and correlation function in this limit are discussed and 
compared to earlier results obtained for systems with homogeneous 
substrate. 

PACS numbers: 68.15.+e; 68.08.Bc 
Keywords: Wetting, Adsorption 
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1 Introduction 



Adsorption phenomena taking place at chemically patterned and geometri- 
cally sculptured substrates have been a topic of growing interest in recent 
years [1-8]. Theoretical research in this field has been stimulated on one side 
by rapid development of experimental methods of imprinting solid substrates 
with structures of sizes reaching nanometer scale [9-12], and on the other by 
attractive future applications, e.g. in the context of micro- and nanofluidics 
or nanotubes [13-14]. 

From theoretical point of view the interfacial morphologies and surface 
phase transitions at inhomogeneous substrates are most often studied via 
methods based on mean-field approach which correctly describe the system's 
behaviour provided intermolecular forces are long-ranged or the system's di- 
mensionality is large enough. On the other hand, fluctuations often mod- 
ify the system's properties, which is best seen in the context of 3d wetting 
with short-ranged forces [TH] or filling transitions in wedge-like geometries 
[2]. Additional interesting features emerge in the case of two-dimensional 
systems with chemically structured substrates, where interfacial fluctuations 
were shown to generate divergences in thermodynamic functions, in partic- 
ular - the point tensions (see [16-17]). Exact results dealing with the case of 
inhomogeneous substrates are still rare though (see [16,18-19]). 

In this paper we are concerned with fluctuation effects accompanying 
adsorption at a one-dimensional planar substrate equipped with a single 
chemical impurity of width 2L. The system under study consists of a two- 
dimensional fluid in a thermodynamic state infinitesimally close to bulk 
liquid-vapor coexistence and exposed to the aforementioned substrate, see 
Fig.l. The chemical structure imposes non-uniformity of the adsorbed liquid- 
like layer, whose morphology and fluctuations are investigated theoretically 
in the present paper. 

The system's temperature is fixed at T = Twi, the critical wetting tem- 
perature of homogeneous substrate type 1, of which the chemical hetero- 
geneity is composed, see Fig.l. The rest of the substrate, which is referred 
to as type 2, remains non-wetted by the liquid; the corresponding wetting 
temperature Tw2 is above Twi- The behaviour of the system in the limit 
of macroscopic inhomogeneity width 2L — > oo is of our main interest here, 
since it corresponds to a typical experimental setup in which the heterogene- 
ity's size is much larger than all correlation lengths present in the system. 
One also notes that the fluctuations' influence on the system is most pro- 
nounced in this limit. On phenomenological grounds, confirmed by earlier 
theoretical predictions [18,20] one may expect that the mean position of the 
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Figure 1: A planar substrate equipped with a chemical inhomogeneity is 
exposed to fluid at bulk liquid-vapour coexistence. A layer of a liquid-like 
phase is adsorbed near the substrate. The interface at y = l(x) separates 
the adsorbed liquid-like layer from the gas phase which is stable in the bulk. 
Type 1 substrate is wetted by the liquid, while substrate type 2 remains 
non-wet. 



interface above the inhomogeneity centre, i.e. at x — 0, as well as interfa- 
cial roughness and parallel interfacial correlation length diverge for L — ► oo. 
On the other hand, as was demonstrated in Ref.|17|. these quantities re- 
main bounded above the chemical impurity's boundaries at x — ±L. It is 
then natural to ask about the behaviour of the characteristic local lengths 
in the intermediate region \x\ e]0,L[, i.e., how strongly the fluctuations are 
suppressed by the non-wet parts of the substrate depending on the distance 
from the heterogeneity's boundary. 

Our approach is based on a continuum interfacial model and allows for 
exact derivation of the equilibrium probability distribution function Pr(h) of 
finding the non-uniform two-phase interface at height lr above a fixed point 
x = T located at the inhomogeneity, as well as the two-point distribution 
function -Pri ,r 2 (^ri ; ^r 2 ) ■ 

The outline of this paper is as follows: in Section 2 we formulate the SOS 
model for the system under study and describe the methods we use to ob- 
tain our results. In Section 3 we derive and discuss the one-point probability 
distribution function -Pr(^r) as well as the two-point probability distribution 
function Pn.r 2 (^ri> ^r 2 )- Section 4 contains calculation of the average inter- 
facial shape (Zr) as well as interfacial roughness (/p) corresponding to the 
limit of macroscopic domain width. In Section 5 the two-point probability 
distribution function is utilized in order to obtain an expression for the in- 
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terfacial correlation function, which depends on the two coordinates I\, T 2 . 
In Section 6 the results are summarized and compared to earlier predictions. 



2 The SOS model 



The present study is based on a continuum interfacial SOS-type Hamiltonian 
model, within which one characterizes the system's state is by a single- valued, 
non-negative function l(x) describing position of the interface above the sub- 
strate at location x. As was demonstrated in Ref.[2T] for the homogeneous 
substrate case, close to wetting temperature this approach captures all the 
system's properties relevant at length scales much larger than bulk correla- 
tion lengths, inter alia, it predicts the same set of surface critical indices as 
those obtained within 2D Ising model [22]. 

For the case of bulk fluid infinitesimally close to liquid- vapour coexistence, 
i.e. for = /ip , /io being the chemical potential at coexistence, the SOS 
Hamiltonian has the following structure 



x 



-X 



a / dl \ 2 



H[l] = dx — I — ) +V(x,l) 



2 \dxJ 



(2.1) 



where a is the interfacial stiffness parameter, and V(x, I) denotes the inter- 
facial potential energy favouring the interface located near the substrate. In 
the present analysis V(x,l) is taken as 

V(x, I) = 0(L - \x\)Vi(l) + e(\x\ - L)V 2 {1) , (2.2) 

where Q(x) is the Heaviside function. In the case of periodic boundary 
conditions l(—X) = /(X), the system's partition function represented by the 
path integral Z(X, L) = J T)le~ n ^ [17,21] can be expressed as 

Z{X,L) = J dl_dl 1 dl 2 Z 2 (l 1 ,L,X-L)Z 1 (l 2 ,l 1 ,2L)Z 2 (l_,l 2 ,X-L) , (2.3) 

where /_ = Z(±X), and li = l(—L), l 2 = l(L) are positions of the interface 
at the inhomogeneity boundaries. The expression for the partition function 
Zi(yi,y 2 , X) corresponding to homogeneous substrate of type i = 1,2, of 
length A, and with interface endpoints fixed at the boundaries at heights yi, 
y 2l respectively, was derived by Burkhardt (see Ref.[2T]) for the case in which 
the potential V^(/) acts only at very short distances. In this paper we are 
concerned exclusively with this case. The partition function Z(X, L) given 
by (|2.3jl was evaluated for arbitrary temperatures in Ref.jEj, where it was 
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then used to discuss the excess point free energy as function of temperature 
and the width L. The probability distribution functions and their moments, 
which can be obtained within this model for the special case T = Ty/i have 
not been discussed so far. 



3 Probability distribution functions 



3.1 One-point probability distribution function 

The one-point probability distribution function of finding the interface at 
height /r above a point Y of the substrate is given by 

Pr(Zr) = (S(l(x = T) - l r )) , (3.1) 

where (....) denotes averaging with the Boltzmann weig ht e~ H V\ with H[l] 
given by Eq.fJSH}; the factor (/c^T) -1 is included into the Hamiltonian. For 
T e] — L,L[ the above formula can be rewritten as 

Pr(Zr) = L ) f dL f dh f dkZ 2 (k,L,X - L) x (3.2) 

xZ 1 (l r , h, L + Y)Z 1 {h, l v , L - T)Z 2 (L, l 2 , X - L) . 

The partition functions Z i {yi ) y 2l \) in Eq. (|3.2j) have the following spectral 
representations in terms of orthonormal sets of solutions {i^n^} to the corre- 
sponding Schrodinger equations with eigenvalues E„ (see Ref.|21|): 

Zifayt, A) = J2^n(yi)^ > (y2)e- E " )x . (3.3) 

n 



For temperatures below T W2 the spectrum {En } corresponding to type 2 
substrate contains a negative lowest eigenvalue Eq . It follows, that replacing 
Z 2 (k,L,X-L) and Z 2 {l_,l 2 ,X-L) in Eq.lQ by ^ ] {h)^* {l^)e~ E o^ x - L ) 
and tpQ (l^)ip^* (l 2 )e~ E a ( X ~ L \ respectively, is equivalent to neglecting terms 

of the order e E o ( X ~ L \ which vanish in the limit of infinite substrate size 
X —y oo. Bearing this in mind one can straightforwardly perform the inte- 
gration over /_ in Eq. (13.21) to obtain 



p(l) = I dh I dk^jh^^Z^lr, h, L + Y)Z 1 {1 2) jr, L — T) 
fdhfdk^ih^i^Z^h^L) ' 
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where the limit X — ► oo has been taken. The remaining integrals can be 

tig the sp 
Twi (see Ref.pT]). namely 



(2) 

performed after inserting the specific forms of and Z 1 corresponding to 
T 



4 2 \y) = V^e Tn 



Z 1 (y 1 ,y 2 ,X) 



a 
2^A 



e -<r(V2~yi? /2A _|_ e -a( yi +y 2 ) 2 /2\ 



(3.5) 



(3.6) 



It is convenient to introduce the 
A+ 



A + = y/(L + T)/£ 2 , 
^r/(2^2±)j where £2 and £ 2 ± are the parallel and 



where r 2 = -\/2oJt[ 2 ~ (T - T W2 ) 
following dimensionless variables: Ao 

A_ = v/(L - r)/6, y 

perpendicular correlation lengths describing interfacial fluctuations above in- 
finite and uniform substrate of type 2. The quantities £ 2) ^2± can be expressed 
in terms of the stiffness parameter a and the eigenvalue [21], and both 
diverge as T approaches T W2 , which is accompanied by the eigenvalue 
increasing to zero. In terms of the above dimensionless parameters the dis- 
tribution Pr(^r) takes the following form 



Pr(y) 



e A °Erfc(A ) 



+ -^A - 2A^e A °Erfc(A 



-1 



(3.7) 



x<jErfc( -J- + X ^ 



+Erfc(^ + A_ 



e- % Erfc( - j- + A_) + Erfc(^- + A_ 
4y Erfc(^- + A_) + Erfc( - j- + A_ 



+ 



where Erfc (a?) stands for the complementary error function 



Erfc(x) 



Erf(x) 



~ t2 dt 



(3.8) 



The probability distribution function Pr(^r) is a decreasing function of Zr at 
any point T, see Fig. 2. Asymptotically this decay is of exponential type with 
characteristic decay length set by £ 2 ±- 



3.2 Two-point probability distribution function 

The two-point probability distribution function -Pri,r 2 (^ri, ^r 2 ) of finding the 
interface at points (r 1; /n) and (r 2 ,/r 2 ) i n the (x,y) plane 

PrinilrM = (Wis = r i) - hMKx = T 2 ) - / r J) (3.9) 
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Figure 2: The probability distribution function Pr(y) plotted for L/^ 2 = 100 
and for two particular choices of T. The solid curve corresponds to the 
inhomogeneity center T = 0, while the dashed one to the close vicinity of the 
inhomogeneity's boundary Y/L = 0.99. Approaching T = is accompanied 
by flattening the distribution function's profile, which indicates growth of the 
interfacial fluctuations' magnitude. 



with Yi G (-L, L) can be written in the form 

Pr 1 ,r.2^r 1 Jr 2 ) = 

j dh j dl 2 4 2 \h)4 2) *(l 2 )Z 1 (l ri ,l 1 ,L + Y 1 )Z 1 (l r2 ,l ri ,Y 2 -Y 1 )Z 1 (l 2 ,l T2 ,L-Y 2 ) 
J dh J dl 2 ^'£\l 1 )^ ) *{h)Z l (l 2 , h, 2L) 

(3.10) 

where we assumed Y 2 > Yi and passed to the limit X — > oo. After inserting 
the specific form of Z\ corresponding to T = Twi, Eq. (j3.6j) into the above 
integrals one obtains 



Pr 1 ,r 2 (yuy2) = e A °Erfc(A ) + -^A - 2Age °Erfc(Ao 



7T 



-1 1 



7T 



Ar 



x (3.11) 



xe A i+ +A H Erfcf - 



m_ 

Ai+ 



+ Ai- 



-2y 



1 + 



Erfcf^ + Ai+U 2?/1 ix 



VA X _ 



. (V2-V1) 



(V2+VlY 



Ar -)-e Ar 



x ^Erfc(-^ + A 2 _)e-^+Erfc(^ + A 2 _)e^^) (t 

where A 1+ = y/ {L + T {) / A 2 _ = y/(L-T 2 )/&, y x = Z ri /(26±), 
Vi = /r 2 /(26±), Ar = (r 2 -ri)/^ 2 - A plot of P Vl ,r 2 (2/1,2/2) for a particular 
choice of system parameters is provided in Fig. 3. 
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Figure 3: The probability distribution function Pr 1 ,r 2 (z/ij 2/2) plotted for 
L/& = 100, Ti/L = 0.1, T 2 /L = 0.2. 



In the next step we shall concentrate on deriving the moments (y) and (y 2 ) 
of the probability distribution function Pr(y) as well as the moment (2/12/2) 
of Pr 1 ,r 2 (l/i ) 2/2) in the asymptotic regime L/£ 2 ^ 1 an d (L — r)/£ 2 ^> 1- 
For this purpose we first note that the arguments of the error functions in 
Pr(y) and -Pri,r 2 (2/i) 2/2) are of the order ^jL/^ 2 (i-e., very large compared to 
1) at any fixed y (or yi, y 2 in the case of Pr u r 2 (yii 2/2))- On the other hand 
for those values of y, yi, y 2 for which the error functions' arguments become 
of the order unity, Pr and -Pri,r 2 are of the order ~ e~ L ^ 2± , see Eq.(3.7), 
(3.11). It follows that the contributions to the integrals determining the 
aforementioned moments that are relevant in the large L regime come only 
from those values of y, y\, y 2 , for which the arguments of the error functions 
are large. One may therefore utilize the asymptotic expansion of Erfc(x) 
around infinity in order to compute the quantities of interest. 
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4 Mean interfacial shape and roughness 



The mean equilibrium interfacial shape (Zr) is given by the first moment of the 
probability distribution function, Eq. (|3.7jl . which we evaluate in the regime 
of macroscopic droplet height, i.e., for L/£ 2 ^ 1 an d for (L — r)/£ 2 3> 1 by 
substituting Erfc(±y/A ± + A±) with ^=>^ e ~ i±y/X±+X±)2 • lt follows that (*r) 
fulfills the following relation 



(lr) I T 2 

for — > oo. The above expression corresponds to the upper root of an 
ellipse, whose long axis of length L is oriented along the substrate, and the 

short axis' length equals y^- Note that the above asymptotic interfacial 

shape does not depend on r 2 . We may now compare (Zr) in Eq. M.ljl to the 
droplet shape calculated for a system consisting of an interface fluctuating in 
presence of a homogeneous 1-dimensional substrate in the case where the in- 
terfacial endpoints are pinned (see Ref . |.2 lj ) . These shapes appear exactly the 
same, which means that the net influence the substrate at x G] — oo, —L[ and 
x e]L, oof exerts on the interfacial morphology above the chemical impurity 
is equivalent to pinning the interface at the inhomogeneity boundaries close 
to the substrate. This conclusion refers to the macroscopic limit L — > oo. 
One also notes, that the scaling interfacial shape in the form of an ellipse 
was obtained for the case of long-ranged van der Waals type interactions via 
mean field approximation, see Ref. [23] ■ From Eq. (|3.7|) one finds the droplet 
height (Zq) divergent as \fh in the considered limit of large L. This may be 
compared to mean-field results [18,20] which predict logarithmic divergence 
of (Z ). As was argued in Ref. [IB], the mean field theory gives correct predic- 
tion regarding (l ) for d=3 system in which interfacial fluctuations are taken 
into account. From Eq. (l4.1j) one concludes, that for T = jL, where 7 G (0, 1), 

(Zr) diverges as y/~L, while in the case T = L — ST, ST being small as 

compared to L, but much larger than £2, one obtains (Z r ) ~ ^(ST) 1 / 2 for 
large L. 

The moment (Zp) can be computed for large L analogously to (Zr). How- 
ever, one may also derive this quantity without reference to any asymptotic 
expansion and for arbitrary values of L. This calculation is sketched in the 
Appendix. It appears, that for L/£ 2 3> 1, (Zp) obeys the following scaling 

(lr) 1 / r 2 \ 
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It follows, that the interfacial roughness (Zp) — (Zp) 2 diverges as -(| — -)(1 — 
^j), which is again the same formula, as one obtains by considering a homo- 
geneous substrate with a fluctuating interface pinned at the boundaries. 



5 Interfacial correlation function 



The SOS model discussed in the present paper gives analytic insight into the 
structure of the correlation function in the limit L/^ 2 1- Let us recall that 
in the present analysis T = T W i which means that type 1 substrate is wetted 
by the liquid. This implies that in the homogeneous case the corresponding 
correlation length £i is infinite. On the other hand the fluctuations are sup- 
pressed by the presence of the outer, semi-infinite non-wet substrate parts. 
Whether for finite L there exists a correlation length ri,r 2 ) describing 
the magnitude of correlations between interfacial fluctuations above points 
Fi and T 2 is not obvious at all. Below we show, that no such quantity arises 
in the regime L/£ 2 ^> 1 and that the correlation function diverges linearly 
for — ► oo provided the points Ti, T 2 are chosen so that (L — Tj) £ 2 
for i = 1,2. 

The two-point interfacial correlation function G(Ti,T 2 , L) is defined by 
G(r 1 ,r 2 ,L) = (Zr 1 Zr 2 )-(Zr 1 )(/r 2 ), (5.1) 

where (lrjr 2 ) = I dlr 1 dlr 2 lr 1 lr 2 Pr 1 ,r 2 (lr 1 Jr 2 )- 

To evaluate (lrjr 2 ) we follow the scheme devised in the preceding sections 
and replace the error functions in Pr 1 ,r 2 ^r 1 Jr 2 ) with their asymptotic at 
infinity. Integration over one of the variables is than straightforward, while 
the second integral is expressed in terms of the hypergeometric function 2 -F\ 
[24] leading to 



/ \ 1 4 1 r A 1+ A 2 _Ar 2 



2(A? + + Ar)(Al_ + Ar) 



(A? + + A|_ +Af)2 2 lj v2' '2' Ar(A? + + Al_ + Af)/J 

M = e A oErfc(A ) + ^=A - 2A 2 e A oErfc(A ) . (5.3) 



where 



An example of the G plot is provided in Fig. 4. To obtain a more transparent 
expression for the correlation function G(Ti, T 2 , L) we focus on the case where 
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Figure 4: The correlation function G plotted as a function of T 2 at fixed 
r x = and L/& = 100 in the range T 2 e (0, L). 



the points Ti, T 2 are separated by a, distance much shorter than L and L — T 2 , 
i.e., we assume Ar <C \f + and Ar <C A 2 _. Using the asymptotic properties 
of the hypergeomertic function 2 -Pi 



,13 x 7rl 11 „, 1 . 



we obtain the asymptotic behaviour of G(Fi, r 2 , L) in the form 



G(r 1; r 2 ,L) iri/ , fA/ r 2 x 1 / v\ I 



a 



(5.5) 



which shows that under the restrictions imposed on the parameters Ti, T 2 
the decay of correlations is in fact linear in Ar, and that for fixed Ti, T 2 the 
correlations grow linearly as a function of L. For Ti = T 2 one recovers the 
local roughness derived in Section 4, Eq.(4.2). The applicability of Eq. (|5.5|) is 
strongly limited by the aforementioned restrictions. On the other hand also 
in formula ()5.2|) one does not observe the presence of terms of type e _Ar ^, 
which typically describe the decay of correlations in non-critical systems. 



6 Summary 



In this paper we analyzed the properties of a two-dimensional model of ad- 
sorption at a substrate equipped with a single chemical impurity. We were 
concerned with the case, when the inhomogeneity is wetted by the liquid 
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phase of the fluid to which the substrate is exposed, while the rest of the 
substrate remains non-wet. The behaviour of mean local interfacial posi- 
tion, interfacial roughness and correlation function was investigated in the 
asymptotic regime of large inhomogeneity's width 2L, for which the system 
becomes critical and is fluctuation-dominated. The type of divergence of 
these quantities in the considered limit was determined depending on the 
distance from the inhomogeneity's boundary. The results were obtained by 
directly calculating the one- and two-point probability distribution functions 
-Pr(Zr) and i^ 1 ,r 2 (^ri> ^r 2 )j an d then investigating the asymptotic properties 
of their moments. Our conclusions are as follows: 

• Both mean interfacial position (Zr) and roughness ((Z r ) — (Ir) 2 ) 1 ^ 2 di- 
verge in the limit L — > oo at any point T which is separated from the 
inhomogeneity's boundary by a distance much larger than the parallel 
correlation length £ 2 corresponding to homogeneous substrate of type 2. 
For T = jL, where 7 G (0, 1), the quantity (Zr) diverges as 

while for T = L — 5T, 5T being small as compared to L, but much larger 
than £ 2 , one obtains (Zr) ~ (5T) 1 / 2 . The mean local height and rough- 
ness close to the chemical structure's boundaries remain bounded even 
for infinite L, as was already noted in Ref . fl7| . The character of the 
divergence of the local roughness is the same as for (Zr). 

• The limiting forms of the mean interfacial shape as well as interfacial 
roughness are independent of the properties of substrate 2 and are the 
same as those obtained for a system composed of an interface fluctu- 
ating above a homogeneous substrate of fixed length 2L provided the 
interface endpoints at — L and L are pinned at a finite distance from 
the substrate. This means that from the point of view of the system's 
large-scale behaviour in the region above the chemical structure, the 
presence of the non-wetted parts of the substrate can be exactly im- 
itated by a system with a homogeneous substrate upon binding the 
interface at — L and L. 

• The obtained profiles of (Z r ) and ((Z r ) — (Zr) 2 ) 1//2 describe ellipses, whose 
x-axis length equals 2L, and y-axis lengths are or the order a/Z and 
are inversely proportional to y/cr. 

• We investigated the structure of the interfacial correlation function 
G(Ti, T 2 , L) in the asymptotic limit of large L and found no indication 
of presence of terms of type e - ^ 2-1 " 1 ^ (r 2 > Ti), which often describe 
decay of fluctuations in non-critical systems. The correlation function's 
amplitude grows linearly with L, and the correlations decay linearly 
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with T 2 — Ti under the restrictions that T 2 — Ti <^ L and r 2 — I\ <C 
L-T 2 . 



7 Appendix 



In this Appendix we discuss the expression for the moment {I 2 ,) of the prob- 
ability distribution function Pr(y) given by Eq. (|3.7|) . The quantity (Z r ) is 
calculated for arbitrary L and T G (— L, L) by noting that Pviy) is invariant 
with respect to the inversion y — > —y. This allows us to write (Zp) in the 
following form 



(y 2 ) 



dyy 2 



Erfc(Ao) H — ^A e 

7T 



\2 
_A 



2A 2 Erfc(A c 



-i 



x 



(A.l) 



e - % Erfc(-^-+A + )Erfc(-^-+A_)+Erfc(-^-+A+)Erfc(^-+A_ 

The above integrals are performed by substituting U (y) = Erfc ^ — + A + j , 

V{(y) = y 2 e - 4y Erfc( - + A_) , V 2 '(y) = y 2 Edc^ + A_) , and integrating 

by parts. This way one arrives at an integral that no longer involves products 
of the error functions. In the next step we insert the integral representation 
fEq. (|3.8|0 of the error function into the obtained formula and integrate over 
y and the remaining variable. This yields the following expression 

-i 



< y 2)= 2 r Erfc(Ao) 2 Aoe -Ag 

'TV L \/7T 



2A^Erfc(A c 



x 



12 



-A 



x 4 \ 2 r 2 



X 



-r 2 +— 

2 16 



(A.2) 
Erfc(Ao) 



from which it follows, that (y 2 ) is quadratic in T at arbitrary L. 

Utilizing the asymptotic expansion Erfc(x) = -j^.e~ x (1 — + 4§i 
+ one obtains the following formula valid for L/£ 2 S> 1 



L 
27 



r 2 

L 2 



r 2 

L 2 



(A.3) 



which up to the leading terms is equivalent to Eq. ()4.2ft . 

One may also explore the asymptotic regime of Eq. (|A.2J) corresponding 
to *C 1. For this case we obtain 



<*r> = 2^ i 1 + 2A o - ^ A o - 8( i A 4 + T 2 + 0{\l 



(A.4) 
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For L/£ 2 = the result for homogeneous substrate type 2 is recovered. One 
observes that the two leading L-dependent corrections are not influenced by 
the value of T. 

The method we used to obtain (/p) for arbitrary L may also be applied 
to compute higher even moments of the distribution Pr(lr)- However, it fails 
in the case of the odd moments, which we were able to obtain only in the 
asymptotic regime L/^ 2 1, as described in the main text. 
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